Order from Disorder in Graphene Quantum Hall Ferromagnet 
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Valley-polarized quantum Hall states in graphene are described by a Heisenberg 0(3) ferromagnet 
model, with the ordering type controlled by the strength and sign of valley anisotropy. A mechanism 
resulting from electron coupling to strain-induced gauge field, giving leading contribution to the 
anisotropy, is described in terms of an effective random magnetic field aligned with the ferromagnet 
z axis. We argue that such random field stabilizes the XY ferromagnet state, which is a coherent 
equal-weight mixture of the K and K' valley states. The implications such as the Berezinskii- 
Kosterlitz-Thouless ordering transition and topological defects with half-integer charge are discussed. 



Gate-controlled graphene monolayer sheets [H host 
an interesting two-dimensional electron system. Re- 
cent studies of transport have uncovered, in particu- 
lar, anomalous Quantum Hall effect resulting from 
Dirac fermion-like behavior of quasiparticles. Most re- 
cently, when magnetic field was increased above about 
20 T, the Landau levels (LL) were found to split with 
the n = ±1 and n = levels forming two and four sub- 
levels, respectively, as illustrated in Fig[T^. The observed 
splittings were attributed to spin and valley degeneracy 
lifted by the Zeeman and exchange interactions. 

The physics of the interaction-induced gapped quan- 
tum Hall state is best understood by analogy with the 
well-studied quantum Hall bilayers realized in double 
quantum well systems 5]. In the latter, the interaction is 
nearly degenerate with respect to rotations of pscudospin 
describing the two wells. As a result, the states with odd 
filling factors are characterized by pseudospin 0(3) order- 
ing, the so-called quantum Hall ferromagnet (QHFM) ^5] . 
The pseudospin z component describes density imbalance 
between the wells, while the x and y components describe 
the inter-well coherence of electron states. Several differ- 
ent phases [1, 0] are possible in QHFM depending on the 
strength of the anisotropic part of Coulomb interaction, 
controlled by well separation. 

In the case of graphene, with all electrons moving in 
a single plane, the valleys K and K' play the role of 
the two wells in the pseudospin representation with the 
lattice constant replacing the inter-well separation. To 
assess the possibility of QHFM ordering, we note that 
the magnetic length at the 10 — 30 T field is much greater 
than the lattice constant. Thus graphene QHFM can 
be associated with the double-well systems with nearly 
perfect pseudospin symmetry of Coulomb interaction [8, 
0, Our estimate, presented below, yields anisotropy 
magnitude of about 10 /xK at i? ~ 30 T, which is very 
small compared to other energy scales in the system. 

Can some other mechanism break pseudospin symme- 
try more efficiently? Coupling to disorder seems an un- 
likely candidate at first glance. However, there is an in- 
teresting effect that received relatively little attention, 
which is strain-induced random gauge field introduced 
by lordanskii and Koshelev [Tlj . To clarify its origin, let 



FIG. 1: a) Graphene Landau level splitting, Ref.0|, at- 
tributed to spin and valley polarization. When the Zeeman 
energy exceeds valley anisotropy, all A'^ = states are spin- 
polarized, with the V = ill states valley-polarized and the 
= state valley-unpolarized. b) The effect of uniform strain 
on electron spectrum, Ref. [ll|] , described by Dirac cones shift 
in opposite directions from the points K and K' . Position- 
dependent strain is described as a random gauge field, Eq.Q. 



us consider the tight-binding model with spatially vary- 
ing hopping amplitudes. Physically, such variation can 

1% ll3| or chemical dis- 



be due to local strain, curvature 
order. With hopping amplitudes ti for three bond orien- 
tations varying independently, we write 
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r*(q) 
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(1) 

where are vectors connecting a lattice site to its nearest 
neighbors, and u and v are wavef unction amplitudes on 
the two non-cquivalcnt sublattices, A and B. The low- 
energy Hamiltonian for the valleys K and K' is obtained 
at q « iqo, where ±qo are two non-equivalent Brillouin 
zone corners: 
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with a 
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f 2 3 '^^1^''^*'*°''''' where the subscript 

+ (— ) corresponds to K{K') valley. Decomposing a± = 
Uy =F iux, we see that the effective vector potential in the 
two valleys is given by ±(ax,ay). Notably, the field a 
is of opposite sign for the two valleys, thus preserving 
time-reversal symmetry (see Fig[T)D). 
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Here we assume that the gauge field has white noise 
correlations with a correlation length ^, 



n A 



(a,(fc)aj(fc)}fe5«i = a^{k) ^ j e '^'^ a^{r)d^r, (3) 

as appropriate for white noise fluctuations of SU. The 
fluctuating effective magnetic fleld can be estimated as 



5h{x) = dxtty 



a/e 



(4) 



whereby the correlator of Fourier harmonics (dhkSh^k) 
behaves as fc^ at fc^ <C 1. 

Recently, strain-induced effective magnetic field was 
employed to explain anomalously small weak localiza- 
tion in graphene [3| • A direct observation of graphene 
ripples [ij] yields typical corrugation length scale ^ of a 
few tens of nanometers. Estimates from the first princi- 
ples lj| gave ~ 0.1 — 1 T, consistent with the observed 
degree of weak localization suppression. 

Valley K-K' asymmetry of QHFM in the presence of 
the gauge field ([2]) translates into a uniaxial random mag- 
netic field, proportional to ([4]) and aligned with the pseu- 
dospin z axis. We shall see that the effect of random 
gauge fleld is subtle: somewhat counterintuitively, weak 
5h induces ordering in the system, acting as an easy plane 
anisotropy which favors the XY state. This behavior can 
be understood by noting that the transverse fluctuations 
in a ferromagnet are softer than the longitudinal fluctua- 
tions, making it beneficial for the spins to be polarized, on 
average, transversely to the field, as illustrated in Figl^l 
This random field-induced ordering maximizes the energy 
gain of the spin system coupled to 5h. 

For magnets with uniaxial random field this behav- 
ior has been established 15, l3| in high space dimension. 
The situation in dimension two is considerably more del- 
icate [itI . [TsI due to competition with the Larkin-Imry- 
Ma fLIMr[l9l.[20t disordered state. We shall see that the 
anisotropy induced by random gauge field is more robust 
than that due to random magnetic field. (This scenario is 
also relevant for the two- valley QH in AlAs system [2l|.) 

The field-induced easy-plane anisotropy completely 
changes thermodynamics, transforming an 0(3) ferro- 
magnet, which does not order in 2d, to the XY model 
which exhibits a Berezinskii-Kosterlitz-Thouless transi- 
tion to an ordered XY state. The transition temperature 
Tbkt is logarithmically renormalized by the out-of-plane 
fluctuations [22.] . Tbkt ^ Jf^i^ilxY /^b), where Ixy is 
the correlation length. For flelds of the order of 30 T, 
with IxY given by Eq. (|15p below, we obtain Tbkt in the 
experimentally accessible range of a few Kelvin. 

The XY-ordered QHFM state hosts fractional ±e/2 
charge excitations, so-called merons 0. Merons are vor- 
tices such that in the vortex core the order parameter 
smoothly rotates out of the xy plane. There are four 
types of merons Q, since a meron can have positive or 
negative vorticity and the order parameter inside the core 
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b) AE = 6n(x)6h(x) = 

FIG. 2: Random field-induced order in a ferromagnet. The 
energy gained from the order parameter tilting opposite to 
the field is maximal when the spins and the field are per- 
pendicular (a), and minimal for the spins parallel to the field 
(b). Uniaxial random field induces ordering in the transverse 
plane. 



can tilt either in +z or —z direction. A pair of merons 
with the same charge and opposite vorticity is topologi- 
cally equivalent to a skyrmion of charge 2(e/2) = e 0]. 

Turning to the discussion of QH effect, the hierarchy of 
the spin- and valley-polarized states is determined by rel- 
ative strength of the Zeeman energy and the randomness- 
induced anisotropy. Our estimate below obtains the 
anisotropy of a few Kelvin at i? 30 T. This is smaller 
than the Zeeman energy in graphene, — gfiBB ~ 50 K 
at S ^ 30 T. Therefore we expect that ^ = state is 
spin-polarized, with both valley states flUed. (This was 
assumed in our previous analysis [2^ of edge states in 

= state.) In contrast, in highly corrugated samples, 
when the anisotropy exceeds the Zeeman energy, an easy- 
plane valley-polarized j/ = state can be favored. 

While the character oi — state is sensitive to the 
anisotropy strength, the i/ = ±1 states (see FigH^) are 
always both spin- and valley-polarized. Below we focus 
on v ^ ±1 states, keeping in mind that for strong ran- 
domnes our discussion also applies to v — state. 

Zeeman-split free Dirac fermion LL are given by 

i?„=sgn(n)|2n|i/25^±A„ iB^{hcleBf'\ (5) 

with n integer and v ~ ?> x lO^cm/s. Each LL is doubly 
valley-degenerate. Random field ([U couples to electron 
orbital motion in the same way as the external field B, 
producing a local change in cyclotron energy and in the 
LL density. While the random field splits the n ^ LL, 
for n = 0, it does not affect the the single-particle energy 
^ and couples to electron dynamics via exchange effects 
only. To estimate this coupling, we note that the field @ 
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leads to valley imbalance in exchange energy per particle: 

Ae^ . /7r\i/2 
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(s) ' 



where k is the dielectric constant of graphene. 

Let us analyze the graphene QHFM energy dependence 
on the gauge field. We consider a fully spin- and valley- 
polarized ly ~ ~1 state, described by a ferromagnetic or- 
der parameter n — [ni, n2, n^) in the K, K' valley space. 
The valley-isotropic exchange interaction gives rise to a 
sigma model, with the gradient term only [6|: 



(7) 



The valley-asymmetric coupling to 5h in Eq. ^ generates 
a Zeeman-like Hamiltonian with a uniaxial random field. 



Ei[n) = / g5h{x)n^{x)d'^x, g = n 
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where n = 1/2-kI\ is the electron density. 

We estimate the energy gain from the order param- 
eter n(x) correlations with the random field, treating 
the anisotropy ([8]) perturbatively in 5h. Decomposing 
n(a;) = n + 5n{x) and taking variation in (5n, we obtain 



(5h 



(z — n(n • z))Sh. 



Substituting the solution for Sn into the energy func- 
tional (O, (IHl), we find an energy gain for n of the form 



6E^-X / {l-ni{x))d^x, A 
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k 
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(9) 



where averaging over spatial fluctuations of 5h is per- 
formed. This anisotropy favors the XY state, = 0. 
Qualitatively (see FiglJ), the fluctuations due to 6n tilt- 
ing towards the z-axis minimize the energy of coupling 
to the uniaxial held when n is transverse to it. 

Now, let us compare the energies of the XY and the 
Larkin-Imry-Ma state 1^, 23| ■ In LIM state the energy 
is lowered by domain formation such that the order pa- 
rameter in each domain is aligned with the average field 
in this domain. Polarization varies smoothly between do- 
mains, and the typical domain size L is determined by 
the balance between domain wall and magnetic field en- 
ergies. In our system, the LIM energy per unit area is 



£lim 



L2 



J 



(10) 



where ^{L) is typical flux value through a region of size 
L. To estimate $(i) we write the magnetic flux through 
a region of size L as an integral of the vector potential 
over the boundary, which gives 



$(L) = j 



at (x) dxi 



(11) 



Minimizing the LIM energy (fTO|) . we find 



Slim 



9^ 



(12) 



Comparison to the XY anisotropy A —g^o^jJ^ gives 
(^)'' ^^-^-1^^ (13) 



J' 



£lim g^a'^ 



where the flux through a region of size ^ where random 
field does not change sign. Interestingly, the ratio (flS]) 
does not depend on the external magnetic field. There- 
fore, at weak randomness, when the random field flux 
through an area is much smaller than the flux quan- 
tum, the ordered XY state has lower energy than the 
disordered LIM state. 

In the opposite limit of strong randomness spins align 
with the local 6h, forming a disordered state. It is in- 
structive to note that for a model with white noise corre- 
lations of magnetic field, rather than of vector potential, 
the ratio (fT3)l is of order one. In this case the competition 
of the LIM and the ordered states is more delicate. 

A different perspective on the random-field-induced or- 
dering is provided by analogy with the classical dynamics 
of a pendulum driven at suspension 2^. The latter, when 
driven at sufficiently high frequency, acquires a steady 
state with the pendulum pointing along the driving force 
axis. As discussed in Ref. [1^, this phenomenon can be 
described by an effective potential UcS obtained by av- 
eraging the kinetic energy over fast oscillations, with the 
minima of UcS on the driving axis and maxima in the 
equatorial plane perpendicular to it. This behavior is ro- 
bust upon replacement of periodic driving by noise [26| . 
Our statistical-mechanical problem differs from the pen- 
dulum problem merely in that the Id time axis is replaced 
by 2d position space, which is inessential for the validity 
of the argument. The resulting effective potential is thus 
identical to that for the pendulum, with the only caveate 
related to the sign change Ucs — C4ff in the effective 
action, as appropriate for transtion from classical to sta- 
tistical mechanics. Thus in our case the minima of Ucs 
are found in the equatorial plane, in agreement with the 
above discussion. 

The easy-plane anisotropy ([9|) can be estimated as 



X/n 



0.1 - 10 K/particle, (14) 



where (5ft, 0.1 — 1 T, ^ 30 nm and _B ^ 30 T was used. 
Since this is smaller than the Zeeman energy, we expect 
that the easy-plane ferromagnet in the valley space is 
realized at v = ±1, while v ~ state is spin polarized 
with both valley states filled. 

The out-of-plane fiuctuations of the order parameter 
are characterized by the correlation length 



XY 



10 £i 



(15) 
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for the above parameter values. The length Ixy sets a 
typical scale for order parameter change in the core of 
vortices (merons) as well as near edges of the sample and 
defects which induce non-zero z-component. 

To measure the correlation length Ixy one may use 
the spatial structure of = wavefunction. Since the 
K(K') electrons reside solely on either A or B sublattice, 
the order parameter z-component is equal to the density 
imbalance between the two sublattices. The latter can 
be directly measured by STM imaging technique. 

Finally, we briefly outline the calculation of QHFM 
valley anisotropy for a pure graphene sheet (the details 
will be published elsewhere). Let us compare the energies 
of state 1, in which only the valley K (or K') Landau level 
is occupied, and state 2, with electrons in an equal-weight 
K-K' superposition state. Since the K {K') electrons at 
1^ = reside on A {B) sublattice, the energies per particle 
in the Hartree-Fock approximation are given by 



reA 



V{r) = -, 



reA,B 



where vq is unit cell volume. (Here we take z = to be a 
site of the A sublattice.) We approximate the energy dif- 
ference El — E2 by the Fourier harmonic of the Hatrtree- 
Fock energy density at the wave vector |Q| 1/a, where 
a ~ 0.14nm is the graphene lattice spacing: 



AE = E, 



E. 



With QIb es/a > 1 at B ~ 30T, the integral yields 
27 



AE: 



512.3 XU)'^^"^°^^ (16) 



B 



indicating that the anisotropy is negligible. 

We note that the situation is completely different 
for higher LL. Goerbig et al. [§| pointed out that the 
Coulomb interaction can backscatter electrons of K and 
K' type at LL with n ^ 0, which leads to a much stronger 
lattice anisotropy of the order q/Ib- This effect is absent 
for the zeroth LL due to the fact that K and K' states 
occupy different sublattices. 

In summary, we studied the valley symmetry break- 
ing of graphene QHFM. We considered the coupling of 
the strain-induced random magnetic field and found that 
it generates an easy-plane anisotropy, which is much 
stronger than the symmetry-breaking terms due to lat- 
tice. The estimates of the field-induced anisotropy sug- 
gest that the random field may be a principal mechanism 
oi K — K' QHFM symmetry breaking. The easy-plane 
ordered state is expected to exhibit BKT transition at 
experimentally accessible temperatures and half-integer 
charge excitations. 
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